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[i]  The  seasonal  variability  of  salinity  transport  in  the  Indian  Ocean  is  investigated  using 
the  high-resolution  global  HYbrid  Coordinate  Ocean  Model  (HYCOM).  Mechanisms 
and  physical  parameters  that  control  the  salinity  budget  are  examined.  Results  show  the 
influence  of  freshwater  forcing  and  zonal  advection  as  the  dominant  mechanisms  of  sea 
surface  salinity  (SSS)  variability.  Precipitation  is  highest  in  the  eastern  Bay  of  Bengal 
(BoB),  where  it  shows  seasonal  variation,  and  in  the  south  equatorial  eastern  Indian  Ocean 
(ElO),  where  it  was  consistently  high  year-round.  These  patterns  result  in  significant 
seasonal  variation  in  the  SSS  in  the  BoB  and  almost  no  variation  in  the  EIO.  Zonal  SSS 
transport  was  higher  than  meridional  SSS  transport  with  the  strongest  seasonality  observed 
along  the  Sri  Lankan  region.  Results  of  depth-integrated  transport  show  northward  salt 
transport  in  the  bottom  layers  and  a  southward  salt  transport  in  the  surface  layers.  The 
4  year  mean  net  flux  of  depth-integrated  salt  transport  was  southward  (-154.8  x  \0t}  kg/s 
to  -552.4  x  I06  kg/s)  at  all  latitudes  except  at  20°N,  where  it  was  northward  (396  *  1 06  kg/s). 

Transport  generally  increases  southward  with  the  highest  transports  occurring  in  the 
south  (1 0°S-35°S)  and  a  maximum  at  30°S.  Analyses  of  meridional  Ekman  volume  and  salt 
transport  show  a  predominantly  southward  transport,  an  indication  of  the  strong  influence 
of  SW  monsoonal  winds.  It  is  anticipated  that  this  study  will  be  useful  in  computing  salt 
transport  using  satellite-derived  salinity  data  from  the  European  Space  Agency  (ESA)  Soil 
Moisture  and  Ocean  Salinity  (SMOS)  and  NASA  Aquarius  salinity  missions. 

Citation:  Nyadjro,  E.  S.,  B.  Subrahmanyam,  and  J.  F.  Shriver  (201 1),  Seasonal  variability  of  salt  transport  during  the  Indian 
Ocean  monsoons,  J.  Geophvs .  Res .,  116,  C08036,  doi;  10. 1029/20 1 1JC006993. 


1.  Introduction 

[:]  The  variability  and  transport  of  salinity  are  important 
in  the  understanding  of  the  global  hydrologic  cycle  and  the 
earth’s  climate,  as  well  as  oceanic  processes  such  as  circu¬ 
lation  and  heat  storage  [Delcroix  et  al. ,  1996;  Wig  and 
Perigaud,  2007].  Except  on  geologic  time  scales  and  the 
consequence  of  climate  change,  net  ocean  salinity  does  not 
change.  The  differences  observed  are  usually  the  result  of 
the  influence  of  precipitation  (P),  evaporation  (E),  river 
runoff  (R),  advection  and  mixing  [ Wijjfels  et  al ,  1992;  Han 
and  McCreary ,  2001;  Prasad  and  Ikeda ,  2002a].  The 
magnitudes  of  E,  P  and  R  define  the  surface  freshwater  flux 
of  a  region  and  they  show  significant  variability  in  the 
Indian  Ocean.  This  is  as  a  result  of  the  seasonal  reversal  of 
winds,  which  causes  monsoon  cycles  to  vary  greatly  thereby 
affecting  the  distribution  and  variability  of  sea  surface 
salinity  (SSS)  [Levi t us,  1988;  Nyadjro  et  al ,  2010].  In  the 
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determination  of  the  conservation  of  salt  in  the  surface 
mixed  layer,  E  and  P  arc  important  parameters  that  are 
incorporated  in  the  mathematical  expressions  of  the  salinity 
budget  [Delcroix  and  Henin ,  1991;  Foltz  and  McPhaden , 
2008].  The  magnitude  of  variability  in  E  and  P  especially, 
and  that  of  other  salinity  tendency  terms  such  as  mixing, 
vertical  and  horizontal  diffusions,  may,  however,  not  lead 
to  an  analogous  magnitude  in  variability  of  salt  transport  as 
the  processes  affecting  the  tendency  terms  vary  in  time  and 
space  [Dour ado  and  Caniaux ,  2003]. 

[3]  Salinity  variability  in  the  Indian  Ocean  can  be  influ¬ 
enced  by  water  circulation  and  ocean -atmosphere  coupled 
interactions  such  as  the  Madden -Julian  Oscillation  (MJO) 
[Joseph  and  Freeland ,  2005].  As  a  consequence  of  the 
differences  in  E-P  and  oceanic  processes  such  as  advection 
and  vertical  mixing,  salinity  varies  among  the  Indian  Ocean 
basins.  Surface  salinity  is  higher  in  the  northwestern  end 
as  a  result  of  E  exceeding  P  in  the  Arabian  Sea  (AS)  and 
lower  in  the  northeastern  end  due  to  higher  precipitation 
in  that  region,  river  runoff  into  the  Bay  of  Bengal  (BoB) 
and  low-salinity  water  from  the  Indonesian  Through  flow 
(1TF)  [Piola  and  Gordon ,  1984;  Donguv  and  Meyers ,  1996; 
Joseph  and  Freeland ,  2005].  For  example,  the  Brahmapu¬ 
tra,  Ganga,  Irrawaddy,  Godavari  and  Mahanadi  rivers, 
which  are  among  the  world’s  50  largest  rivers,  flow  into  the 
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BoB  [Sengupta  et  al. ,  2006].  It  is  estimated  that  4700  and 
3000  kmVyr  of  freshwater  enters  BoB  via  local  precipita¬ 
tion  and  river  discharge  respectively;  while  3600  kmVyr  is 
lost  via  evaporation  [Rajamani  et  al .,  2006].  On  the  whole, 
the  Indian  Oecan  is  considered  an  evaporative  basin  as  it 
loses  about  1.392  *  104  km3/yr  of  freshwater  in  net  evap¬ 
oration  to  the  atmosphere  [Wijffels  et  al. ,  1992].  The  vari¬ 
ation  in  the  flux  parameters  therefore  ean  set  up  seasonal 
differences  in  the  initiation,  sustenance  and  variability  of 
salinity  transport  in  the  Indian  Ocean  especially  given  the 
dynamics  of  processes  associated  with  the  monsoons. 

[4]  Palcoceanographie  records  have  shown  that  the 
strength  of  monsoons  varies  on  several  time  scales  and  their 
effects  on  oceanographic  conditions  and  processes  could  as 
well  be  varied  [ Webster  et  al .,  1998;  Murtugudde  et  al ., 
2007].  Nevertheless,  the  causes  and  patterns  in  salinity 
distribution  and  variability  are  quite  consistent.  For  exam¬ 
ple,  although  one  of  the  main  indicators  of  the  strength  of  a 
monsoon  is  the  magnitude  of  precipitation  (which  signifi¬ 
cantly  affects  salinity  values),  Vinayachandran  [2004] 
contends  that  the  strength  and  duration  of  monsoonal 
winds  were  more  important  for  AS  cooling  (and  by  infer¬ 
ence  salinity  variation)  than  the  amount  of  precipitation  into 
the  BoB  and  subsequently  into  the  AS. 

[5]  In  the  Indian  Ocean,  salinity  contributes  to  the  for¬ 
mation  of  the  barrier  layer  as  a  result  of  stratification  w  ithin 
the  isothermal  layer,  controls  vertical  mixing,  mixed  layer 
and  warm  pool  dynamics  [Han  et  al. ,  2001;  Durand  et  al ., 
2007;  Hareesh  Kumar  et  al .,  2009].  The  transport  of  salt 
therefore  has  the  potential  of  affecting  these  processes. 
Additionally,  it  is  important  in  the  understanding  of  the 
global  water  budget  which  is  one  of  the  main  forcing 
mechanisms  of  the  thcrmohalinc  circulation.  Increase  in 
surface  salinity  when  aided  by  deep  convection  could  lead 
to  the  formation  of  deep  water  masses.  This  phenomenon 
has  the  potential  of  inducing  a  meridional  density  gradient 
that  affects  meridional  circulation  and  aids  the  transport 
of  salt  [Sevellec  et  al .,  2008;  Czaja ,  2009].  Deep  convection 
could  be  inhibited  by  the  freshening  of  the  surface  waters. 
To  sustain  deep  convection  and  meridional  circulation, 
ocean  currents  transport  salt  to  balance  out  the  freshening 
[Oka  and  Hasumi ,  2006]. 

[6]  Surface  circulation  in  the  Indian  Ocean  undergoes 
seasonal  variability  as  a  result  of  reversal  of  winds,  the 
effects  of  which  is  well  noticeable  in  the  northern  part  of 
the  ocean.  The  Southwest  (SW)  monsoon  occurs  between 
June  and  September  while  the  Northeast  monsoon  (NE) 
occurs  between  November  and  February  .  The  onset  of  the 
SW  monsoon  is  usually  characterized  by  weak  alongshore 
winds  off  the  coast  of  Somalia.  It  then  matures  in  June  and 
peaks  in  the  later  parts  of  July  [Webster  et  al. ,  1999].  During 
the  SW  monsoon,  the  highly  saline  Arabian  Sea  Water  mass 
(ASW),  w  hieh  oecurs  at  a  depth  of  about  0-100  m,  spreads 
with  the  eastward  flowing  SW  Monsoon  Current  (SMC) 
carrying  salty  water  into  the  BoB  [ Morrison ,  1997;  Schott 
and  McCreary ,  2001].  The  Somali  Current  also  flows 
northward  into  the  AS.  During  this  season  too,  the  humid 
southwesterly  winds  from  the  oceans  causes  increased  sur¬ 
face  cooling  and  significant  upwelling  along  the  coasts  of 
the  Arabian  Seas  [Murtugudde  et  al. ,  2007].  This  is  capable 
of  bringing  high  saline  waters  to  the  sea  surface  and  coupled 
with  high  evaporation,  makes  SSS  in  this  region  usually 


higher  than  other  times  of  the  year  [Prasad  and  lkeda , 
2002b].  The  NE  monsoon  is  characterized  by  dry  and 
weak  northwesterly  winds  from  Asian  landmass.  In  response 
to  the  reversed  winds,  the  Somali  current  also  reverses  and 
flows  southwards.  During  this  season,  less  saline  waters  from 
the  BoB  is  earned  by  the  southward  flowing  East  India 
Coastal  Current  (E1CC)  into  the  West  India  Coastal  Current 
(W1CC)  [Vinayachandran  et  al .,  1999;  Prasanna  Kumar 
et  al .,  2004].  The  westward  flowing  Northeast  Monsoon 
Current  (NMC)  also  brings  less  saline  water  into  the  WICC 
which  then  flows  northward  along  the  west  coast  of  India 
into  the  AS  [Schott  and  McCreary ,  2001;  Subrahmanyam 
et  al.,  2011].  This,  however,  does  not  result  in  drastie  drop 
in  AS  SSS  as  the  passage  of  the  dry  NE  winds  results  in 
evaporative  cooling,  which  ensures  a  high  salinity  [Prasad 
and  lkeda ,  2002a]. 

[7]  The  water  circulation  in  the  southern  Indian  Ocean  is 
predominantly  a  eonsequenee  of  a  wind  stress  eurl  forced 
subtropical  gyre  that  occurs  in  the  upper  1000  m  and  made 
up  of  the  Agulhas  Current,  East  Madagascar  Current 
(EMC),  Mozambique  Current  (MC)  and  the  South  Equato¬ 
rial  Current  (SEC).  The  net  flow  of  these  is  southward  and 
dominated  by  the  Agulhas  Current  [Gordon  et  al .,  1992]. 
The  flows  in  the  deep  and  abyssal  sections  of  the  region  are 
dominated  by  the  presence  of  bottom  currents  that  transport 
the  Circumpolar  Deep  Water  (CDW,  occurring  at  ~2000 
3000  m)  and  the  Antarctic  Bottom  Water  (AABW,  occur¬ 
ring  below  ~3800  m)  northward  [ Fu ,  1986].  The  CDW 
enters  the  southwestern  end  of  the  Indian  Ocean  through  the 
Agulhas  Basin  then  into  the  Mozambique  Basin,  its  north¬ 
ernmost  extent  is,  however,  believed  to  be  Davie  Ridge 
(~2.5  km  high)  in  the  Mozambique  Channel  [Toole  and 
Warren ,  1993;  Mantyla  and  Reid ,  1995].  Later  studies  by 
van  Aken  et  al.  [2004],  however,  suggest  that  about  2  Sv  of 
the  top  half  of  the  NADW  flow'  across  the  Ridge  into  the 
Somali  Basin.  This  water  travels  further  north  into  the 
Arabian  Basin,  where  it  gradually  upwclls  into  and  becomes 
part  of  the  overlying  Indian  Deep  Water  (IDW).  The  flow  of 
this  water  is  predominantly  northward  and  its  influence 
causes  northward  salt  and  volume  fluxes.  This  water  which 
typically  oecurs  between  1500  to  3800  m  also  accounts  for 
the  northw  ard  flux  of  w  aters  in  the  intermediate -deep  layers. 
In  the  Arabian  Basin,  high-saline  Red  Sea  and  Persian  Gulf 
waters  vertically  mix  with  the  underlying  waters  thereby 
increasing  the  salinity  of  the  deeper  waters  [Piola  and 
Gordon ,  1984].  The  Red  Seawater  has  been  reported  to 
have  been  seen  as  far  south  as  the  southern  end  of  the 
Mozambique  Channel  [ Grundlingh ,  1985]. 

[#]  Due  to  the  limitation  of  salinity  observations,  basin 
scale  studies  of  salinity  transport  have  been  a  challenge  in 
the  Indian  Ocean.  Many  such  studies  have  relied  heavily  on 
numerical  simulations,  which  sometimes  suffer  especially 
from  the  estimation  of  freshwater  flux  terms  [Murtugudde 
and  Busalacchi ,  1998].  This  has  been  explained  as  arising 
from  the  fact  that  the  spatial  distribution  of  atmospheric 
input  functions  is  non-Gaussian  [ Tomczak ,  1995;  Dourado 
and  Caniaux ,  2003].  Other  challenges  have  been  the  poor 
representation  of  coastal  currents  and  the  dominance  of 
Ekman  drift  over  geostrophy  in  the  western  coast  of  India 
that  leads  to  differences  between  observed  and  modeled  SSS 
[Prasad  and  lkeda,  2002b].  Using  passive  tracers  to  map  the 
exchange  of  salt  in  a  multilayer  upper  ocean  model,  Jensen 
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[2001]  attributed  the  southward  advection  of  low-saline 
water  from  the  BoB  to  Ekman  drift.  His  model,  however, 
did  not  find  a  significant  transport  of  water  from  the  BoB 
into  the  AS  although  Esenkov  et  al.  [2003]  did.  Esenkov 
et  al.' s  [2003]  Miami  Isopycnic  Coordinate  Ocean  Model 
(M1COM)  could  not  capture  the  influence  of  the  pathway 
between  Sri  Lanka  and  India,  known  as  the  Palk  Strait,  in 
the  exchange  of  salt  between  AS  and  BoB  as  this  was  closed 
in  the  model.  In  our  HYCOM  simulations,  the  Palk  Strait 
was  opened  and  allows  water  exchange  between  the  BoB 
and  the  Gulf  of  Mannar. 

[9]  The  bulk  of  the  Indian  Ocean  lies  in  the  Southern 
hemisphere  for  which  data  is  very  sparse.  In  spite  of  its  role 
in  global  climate  phenomena,  the  Indian  Ocean  is  poorly 
sampled  and  not  many  salt  transport  estimates  exist  for  it. 
The  few  studies  available  are  limited  to  a  small  number  of 
sections  (typically  one  or  two  latitudinal  bands)  thereby 
making  the  basin-wide  understanding  of  processes  quite  a 
challenge.  To  the  best  of  our  knowledge  therefore,  this  is 
the  first  basin-scale  depth -integrated  salt  transport  study  for 
the  Indian  Ocean.  It  hopes  to  build  on  earlier  works  by  using 
a  1/12°  high -resolution  HYbrid  Coordinate  Ocean  Model 
(HYCOM)  simulation  during  2003-2006  to  study  salt 
transport  in  the  Indian  Ocean.  It  describes  the  seasonal 
variability  of  salinity  and  the  relative  importance  of  flux 
terms  in  the  advection  process.  It  also  investigates  the 
pathways  of  salt  exchange  as  a  function  of  depth  in  the 
Indian  Ocean.  The  complex  topography  of  the  Indian  Ocean 
offers  bathymetric  constraints  for  which  the  approach  used 
in  the  study  does  not  clearly  highlight  the  role  of  depth 
features  such  as  sills  in  the  transport;  nevertheless,  a  general 
pattern  of  the  salt  transport  is  achieved. 

2,  Model,  Data,  and  Methods 

2.1.  HYCOM  Simulations 

[10]  The  model  used  in  this  study  is  global  HYCOM, 
which  was  developed  from  the  Miami  Isopycnic  Coordinate 
Ocean  Model  (Ml COM)  using  the  theoretical  foundation  set 
forth  by  Bleck  and  Boudra  [1981],  Bleck  and  Benjamin 
[1993]  and  Bleck  [2002].  This  model  has  1/12°  horizontal 
resolution  (~7  km  at  midlatitudes  and  3.5  km  at  the  North 
Pole)  and  32  hybrid  layers  in  the  vertical,  is  isopycnal  in 
the  open  stratified  ocean  and  makes  a  smooth  transition  to 
a  terrain-following  (rr)  coordinate  in  coastal  waters  using 
the  layered  continuity  equation  [Bleck.  2002].  The  model’s 
use  of  varying  vertical  coordinates  allows  it  to  effectively 
simulate  the  characteristics  of  ocean  circulation  as  well 
as  interior  water  mass  distribution.  These  features  are  very 
important  especially  in  the  study  of  transport  of  ocean 
properties  such  as  salinity  as  the  dynamics  arc  strongly 
influenced  by  localized  processes  hence  the  use  of  an 
entirely  single  vertical  coordinate  such  as  depth  will  be 
inefficient  [Chassignet  et  al .,  1996].  The  model  is  config¬ 
ured  on  a  Mercator  grid  from  78°S  to  47°N,  while  north 
of  this  latitude  an  Arctic  dipole  patch  is  used  to  avoid 
the  singularity  at  the  pole.  HYCOM  has  the  capability  of 
selecting  among  several  different  vertical  mixing  submodels 
(e.g.,  K- Pro  file  Parameterization  [Ixirge  et  al .,  1994]  and 
Dynamical  Instability  Model  [Price  et  a/.,  1986])  for  either 
the  surface  mixed  layer  or  the  weak  interior  diapycnal 
mixing.  This  feature  is  essentially  important  for  our  study  as 


salinity  distribution,  especially  in  the  surface  ocean  is  partly 
dependent  on  vertical  mixing. 

[11]  Our  version  of  HYCOM  was  initialized  using 
monthly  mean  temperature  and  salinity  from  the  1/4° 
Generalized  Digital  Environmental  Model  (GDEM3)  cli¬ 
matology  l NAVOCEANO ,  2003]  in  August.  This  experi¬ 
ment  was  then  run  for  13  years  using  climatological 
monthly  mean  wind  and  thermal  forcing  constructed  from 
the  1.125°  European  Centre  for  Medium-Range  Weather 
Forecasts  (ECMWF)  Re-Analysis  (ERA  15)  over  the  1 979— 
1993  time  frame.  In  addition,  6-hourly  variability  from  the 
ECMWF  operational  model  over  the  period  September 
1 994-September  1995  was  added  to  the  climatological 
wind  forcing  (but  not  the  thermal  forcing)  to  add  the  higher- 
frequency  variability  needed  for  realistie  simulation  of  the 
surface  mixed  layer.  In  order  to  keep  the  evaporation  minus 
precipitation  budget  on  track,  the  model  weakly  relaxes  to 
monthly  mean  SSS  from  the  Polar  Science  Center  Hydro- 
graphic  Climatology  (PHC)  SSS.  This  climatology  is  cho¬ 
sen  for  its  accuracy  in  the  Arctie  region  [Steele  et  al .,  2001]. 
The  actual  SSS  relaxation  e-folding  time  depends  on  the 
mixed  layer  depth  (MLD)  and  is  (30  days  *  30  m/MLD  m) 
days,  i.e.,  it  is  more  rapid  when  the  MLD  is  shallow  and 
less  so  when  it  is  deep.  The  SSS  relaxation  is  sufficient  to 
avoid  long-term  drift  in  SSS,  but  not  so  strong  as  to  inhibit 
SSS  anomaly  formation.  The  SSS  relaxation  is  in  addition 
to  the  evaporation -precipitation  budget. 

[12]  Once  the  model  was  determined  to  be  equilibrated 
after  being  spun -up  with  climatological  forcing,  the  simu¬ 
lation  was  then  continued  using  3  hourly  Navy  Operational 
Global  Atmospheric  Prediction  System  (NOGAPS)  winds, 
heat  fluxes  and  precipitation  fields  for  the  period  2003-2006 
as  the  forcing.  This  version  of  the  model  uses  the  NASA 
Goddard  Institute  for  Space  Studies  level  2  (G1SS)  mixed 
layer  scheme  and  includes  monthly  river  runoff  from  986 
global  rivers.  There  is  no  assimilation  of  any  ocean  data, 
including  SST,  and  no  relaxation  to  any  other  data  except 
SSS  to  keep  the  evaporation  minus  precipitation  balance  on 
track. 

[13]  Because  we  are  using  the  global  version  of  HYCOM, 
there  is  no  sponge  layer  across  the  lateral  boundaries  and  the 
model  allows  the  ITF  to  flow  into  the  Indian  Ocean  and  salt 
would  be  conserved.  Also  salt  transport  is  estimated  in 
various  sigma  t  layers  wherein  it  is  expected  that  salt  and 
mass  arc  conserved  and  no  attempt  is  placed  to  estimate  the 
diapycnal  mixing. 

2.2.  Data 

[14]  Climatological  data  sets  on  Evaporation  and  Precip¬ 
itation  for  surface  salinity  variability  analyses  were  obtained 
from  the  Hamburg  Ocean  Atmosphere  Parameters  and 
Fluxes  from  Satellite  data  set  (HOAPS),  Meteorological 
Institute  of  the  University  of  Hamburg,  Germany  [Anderson 
et  al .,  2007].  These  parameters  are  generated  from  data 
obtained  from  the  Special  Sensor  Microwave/lmager  (SSM/ 
1)  radiometer  aboard  the  Defense  Meteorological  Satellite 
Program  (DMSP)  satellites.  The  surface  wind  stress  r  = 
PJCdW  was  computed  from  the  Florida  State  University’s 
(FSU)  Center  for  Ocean-Atmosphere  Prediction  Studies 
(COAPS)  1 0  x  1 0  gridded  pseudostress  (^)  data  with  pa  = 
1.2041  kgW  and  CD  =  0.0015  [Bourassa  et  al .,  2005]. 
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Table  I.  HYCOM  Layers  With  Approximate  Depth 


Layer 

Mean  Thickness  (m) 

Approximate  Depth  (m) 

1 

5 

5 

2 

5 

10 

3 

10 

20 

4 

10 

30 

5 

20 

50 

6 

25 

75 

7 

25 

100 

8 

50 

150 

9 

50 

200 

10 

100 

300 

11 

100 

400 

12 

100 

500 

13 

100 

600 

14 

100 

700 

15 

100 

800 

16 

100 

900 

17 

200 

1100 

18 

200 

1300 

19 

200 

1500 

20 

200 

1700 

21 

250 

1950 

22 

250 

2200 

23 

300 

2500 

24 

300 

2800 

25 

300 

3100 

26 

500 

3600 

27 

500 

4100 

28 

500 

4600 

29 

750 

5350 

30 

750 

6100 

31 

500 

6600 

32 

500 

7100 

2.3.  Methods 

[15]  Hastenrath  [1984]  summarizes  the  approach  of  esti¬ 
mating  meridional  transport  as  (a)  evaluating  the  surface 
salinity  budget  and  subsequently  integrating  with  latitude  to 
obtain  transports  across  latitude  circles,  (b)  directly  evalu¬ 
ating  transports  from  zonal  hydrographic  sections  and  (c) 
estimating  the  divergence  of  zonal  transport  by  a  major 
current  (c.g.,  ACC)  from  meridional  hydrographic  sections. 
We  chose  the  first  approach  as  it  will  enable  us  to  investigate 
the  direct  impact  of  seasonal  reversal  of  winds  and  currents 
on  surface  ocean  salinity  variability  and  transport  and  then 
link  it  up  to  the  subsurface  transport.  Additionally,  the  esti¬ 
mation  of  sea  surface  salinity  transport  from  HYCOM  will 
allow  for  future  comparison  with  that  estimated  using  sat¬ 
ellite-derived  surface  salinity  measurements  from  the  Euro¬ 
pean  Space  Agency  (ESA)  Soil  Moisture  and  Ocean  Salinity 
Mission  (SMOS)  and  the  future  NASA  Aquarius  Salinity 
Mission. 

[16]  SSS  transport  in  the  Indian  Ocean  was  computed 
from: 

Fsvs  =  pvSA  (I) 


v  is  replaced  with  w,  the  zonal  velocity  in  m/s.  The  surface 
salinity  transport  is  then  connected  to  the  subsurface  ocean 
by  computing  depth-integrated  salt  transport  across  latitude 
circles  using  the  HYCOM  salinity  along  a  section.  In  this 
regard,  we  compute  salt  transports  along  5°  latitudinal 
transects  in  the  Indian  Ocean  using  monthly  mean  HYCOM 
v  ariables  and  equations  from  Wijffels  et  al  [1992],  Saunders 
and  Thompson  [1993]  and  Ganachaud  et  al  [2000]: 

z„  L 

F,=  J  J p.v(x,z)S{x,z)dxdz  (2) 

-//  O 

where  Fs  is  the  zonally  integrated  meridional  salt  transport 
in  106  kg/s,  v  is  the  meridional  oceanic  velocity  in  m/s  and 
is  perpendicular  to  the  transect  and  positive  for  northward 
velocity,  S  is  the  spatially  varying  oceanic  salinity,  p  is  the 
density  in  kgW  estimated  from  the  potential  density  of 
the  respective  32  isopycnal  layers  in  the  HYCOM,  x  and  z 
are  the  zonal  and  vertical  coordinates  respectively,  z(i  is 
the  ocean  surface,  H  is  the  water  depth  in  meters  at  a  fixed 
position,  jr  and  L  is  the  zonal  expanse  in  meters  of  the  ocean 
section  at  depth  //.  To  compute  the  salt  transport  as  per 
equation  (2),  information  on  the  vertical  distribution  of 
salinity  is  required.  To  do  this,  we  used  the  interface  depth  of 
the  isopycnal  layers  in  the  HYCOM  simulations,  obtaining 
the  approximate  depths  of  the  layers  as  shown  in  Table  1. 
Using  isopycnal  layers  in  the  open  ocean  allows  transport  to 
follow  surfaces  of  constant  potential  density  in  the  Indian 
Ocean.  It  w  ill  also  allow  the  interpretation  of  transport  results 
w  ith  respect  to  water  masses.  Additionally,  volume  transport 
was  computed  to  aid  the  interpretation  and  analyses  of  salt 
transport  results;  thus 

L 

Fv  =  j  J  v(x,z)dxdz  (3) 

//  o 

where  Fv  is  volume  transport  in  Sv  (1  Sv  =  106  m1  s  ').  For 
estimates  of  transports  via  the  Indonesian  Throughflow 
(1TF),  the  meridional  velocity  v  in  the  above  equations  2  and 
3  is  replaced  with  w,  the  zonal  velocity  and  is  positive  for 
eastward  velocity.  Results  of  the  computations  arc  shown  for 
25°N  to  35°S,  where  the  Indian  Ocean  has  an  east  and  west 
landmass  boundary'. 

[n]  In  an  earlier  study  [Nyadjro  et  al. ,  2010),  Ekman 
component  was  reported  to  play  a  significant  role  in  the 
distribution  of  salt  in  the  sea  surface.  In  this  study,  we 
analyze  the  effect  of  w  inds  in  salt  transport  in  the  Ekman 
layer  by  computing  depth  integrated  meridional  Ekman 
transports  using  COAPS  winds  data.  The  annual  mean 
meridional  Ekman  volume  transport  from  HYCOM  was 
computed  following  similar  methods  of  Levitus  [1988]: 


where  Fsss  is  the  surface  salinity  transport  in  106  kg/s,  p 
is  the  density  of  sea  surface,  1023  kg/m\  S  is  sea 
surface  salinity  (defined  as  *10  ,  hence  the  transport  units 
of  106  kg/s),  v  is  the  meridional  velocity  in  m/s,  A  is  the 
cross-sectional  area  computed  from  h  *  /,  where  h  ~  5  m 
(the  depth  of  the  first  model  layer)  and  /,  the  zonal  extent  of 
the  cross-sectional  area.  For  zonal  surface  salinity  transport, 


(4) 


where  Q 0  is  the  total  horizontal  Ekman  transport,  T\  is  the 
zonal  component  of  the  wind  stress,  /  is  the  Coriolis 
parameter  and  p  is  the  density  of  seawater  from  HYCOM. 
Ve  is  the  meridional  Ekman  velocity  component  and  ~h  is 
the  depth  of  the  Ekman  boundary  layer,  which  wc  chose  to 
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Figure  1.  Mean  annual  (a)  HYCOM  SSS,  (b)  E  (mm/day),  (c)  P  (mm  day),  and  (d)  E-P  (mm/day)  with 
boxes  denoting  the  averaging  subregions  for  surface  salinity  variability'  terms  and  surface  salinity  trans¬ 
port  studies:  Arabian  Sea  (AS),  Bay  of  Bengal  (BoEJ),  Sri  Lanka  (SL),  and  Equatorial  (EQ)  regions. 


be  100  m  (~layer  7  of  HYCOM).  Meridional  Ekman  salt 
transport  was  computed  from  Chereskin  el  al.  [2002]: 

r„  l 

Sf=  j  J  p.V£Sdxd. z  (5) 

//  0 

where  Sg  is  the  meridional  Ekman  salt  transport,  S  is 
HYCOM  salinity. 

[is]  We  compare  our  depth- integrated  salt  transport  from 
HYCOM  with  that  using  the  Simple  Oeean  Data  Assimi¬ 
lation  (SODA)  model  version  2.2.4  reanalysis  output  first 
described  by  Giese  and  Ray  [201 1].  This  version  of  SODA 
combines  an  ocean  model  based  on  the  Parallel  Ocean 
Program  (POP)  version  2.0.1  numerics  [Smith  et  al.,  1992] 
with  hydrographic  temperature  and  salinity  data.  SODA  is  a 
data-assimilated  product  and  has  a  horizontal  resolution  that 
is  on  average  0.4°  *  0.25°  and  with  40  levels  in  the  vertical. 
The  atmospheric  reanalysis  includes  only  surface  observa¬ 
tions  of  synoptic  pressure  and  monthly  SST  and  sea  ice 
distribution  from  the  Hadley  Center  lee  Sea  Surface  Tem¬ 
perature  (HadlSST  v.  1 . 1 )  data  set  [Rayner  el  al .,  2003].  The 
temperature  and  salinity  profile  data  used  in  SODA  was 
obtained  from  the  reecnt  release  of  the  World  Oeean  Data¬ 
base  2009  (WOD09)  [Boyer  el  al .,  2009]. 

3.  Results  and  Discussion 

3.1.  Sea  Surface  Structures  and  Standard  Deviation 

[19]  We  present  E,  P,  E-P  results  for  spatially  averaged 
data  for  the  Arabian  Sea  (AS),  Bay  of  Bengal  (BoB),  Sri 
Lanka  (SL)  and  the  Equatorial  (EQ)  regions  which  respec¬ 
tively  arc  represented  by  the  following  boxes:  AS  (10°N  to 


20°N  and  60°E  to  70°E),  BoB  (12°N  to  18°N  and  83°E 
to  93°E),  SL  (2°N  to  10°N  and  75°E  to  85°E)  and  equatorial 
(5°S  to  15°S  and  70°E  to  90°E)  (Figure  la).  These  study 
regions  are  selected  to  reflect  the  different  climatological 
and  physical  forcings  in  the  Indian  Ocean.  Annual  mean  for 
SSS,  E,  P  and  E-P  are  respectively  shown  in  Figures  1  a—  1  d. 
SSS  is  highest  in  the  AS  and  lowest  in  the  BoB  (Figure  la) 
with  the  difference  mainly  as  a  result  of  E-P  being  greater  in 
the  former  than  the  latter  (Figure  Id).  Although  the  annual 
mean  evaporation  seem  to  be  similar  for  the  two  regions 
(Figure  lb),  precipitation  is  greater  in  the  BoB,  especially 
during  the  SW  monsoon.  Additionally  the  high  river  runoff 
into  the  Bay  also  contributes  to  its  lower  salinity  than  the 
AS  [Sengupta  el  al .,  2006].  The  region  south  of  the  equator 
(10°S-30°S)  shows  the  highest  evaporation  zone  in  the 
Indian  Oeean  (Figure  lb).  Although  the  E-P  in  this  region  is 
similar  to  that  of  the  AS  (Figure  Id),  SSS  is  lower  here 
beeause  of  fresher  water  from  the  1TF  [Sengupta  el  al ., 
2006]  while  the  AS  additionally  receives  high  saline 
waters  from  the  Red  and  Persian  Seas  at  subsurface  depths. 

[20]  The  temporal  variability  of  the  surface  salinity  vari¬ 
ability  terms,  E,  P  and  E-P  arc  shown  in  Figure  2.  On  the 
average,  precipitation  is  highest  in  the  BoB  and  the  EQ 
(Figures  2a  and  2b)  especially  in  the  eastern  ends  of  these 
regions  (Figure  le).  Whereas  in  the  EQ  region  it  is  consis¬ 
tently  high  throughout  the  years,  that  of  the  BoB  shows 
significant  seasonal  variation.  The  annual  mean  and  standard 
deviation  respectively  are  5.59  mm/day  and  5.46  mm/day  for 
the  BoB  and  6.50  mm/day  and  1.87  mm/day  for  the  EQ 
region.  During  the  NE  monsoon  BoE*  precipitation  is  as  low 
as  <  I  mm/day  but  rises  to  about  1 5  mm/day  at  the  peak  of 
the  SW  monsoon  with  a  month  lag  behind  the  precipitation 
peak  in  the  AS.  In  the  AS,  precipitation  shows  seasonal 


5  of  14 


C '08036 


NYADJRO  ET  AL.:  SALT  TRANSPORT  IN  THE  INDIAN  OCEAN 


(08036 


Figure  2.  Temporal  variability  of  surface  salinity  variability  terms:  (a,  b)  precipitation  (P),  (e,  d)  evap¬ 
oration  (E),  and  (e,  0  E-P.  Figures  2a,  2c,  and  2e  are  for  AS  and  BoB  regions,  while  Figures  2b,  2d,  and 
2f  are  for  SL  and  EQ  regions. 


variation  with  a  near  zero  value  throughout  the  NE  monsoon 
season  but  peaks  to  about  5  mm/day  during  the  SW  monsoon 
season  (Figure  2a).  Annual  mean  and  standard  deviation  for 
this  region  are  1 .49  mm/day  and  1 .99  mm/day  respectively. 
Precipitation  in  the  SL  region  (Figure  2b)  depicts  a  monthly 
sinusoidal  pattern  without  a  clear  seasonal  variation  with 
annual  mean  and  standard  deviation  being  3.94  mm/day  and 
2.12  mm/day. 

[21]  The  evaporation  term  shows  seasonal  variation  in  all 
the  study  regions  (Figures  2e  and  2d).  During  the  NE 
monsoon  and  at  the  start  of  the  SW  monsoon,  evaporation  in 
the  four  regions  are  all  in  phase.  At  the  peak  of  the  SWr 
monsoon,  however,  that  of  SL  and  EQ  increases  while  it 
decreases  in  the  AS  and  BoB  regions.  The  annual  mean 
and  standard  deviation  respectively  are  4.68  mm/day  and 
0.67  mm/day  for  the  SL  and  5.11  mm/day  and  0.82  mm/day 
for  the  EQ  region.  The  AS  has  an  annual  mean  of  4.96  mm/ 
day  and  a  standard  deviation  of  1.06  mm/day  while  for  the 
BoB  it  is  4.51  mm/day  and  0.82  mm/day  respectively.  The 
AS  region  (Figure  2c)  shows  annual  bimodal  peaks  occur¬ 
ring  in  June  and  December  with  the  former  coinciding  with 
the  minima  E-P  and  the  latter  with  the  maxima  E-P  in  the 
AS  region  (Figure  2c). 

[22]  The  largest  variation  in  the  E-P  term  is  seen  in  the 
BoB  (Figure  2e)  where  the  temporal  distribution  generally 
is  a  mirror  image  of  the  P  variations  in  the  BoB,  an  indi¬ 
cation  of  the  controlling  role  of  precipitation  in  this 
region  and  the  effect  it  has  on  salinity  flux.  The  annual 
mean  E-P  is  -1.07  mm/day  with  a  high  standard  deviation 
of  5.79  mm/day.  The  net  flux  ean  reaeh  -15  mm/day  in 
July,  implying  an  input  of  freshwater  into  the  basin.  Along 
the  EQ  region  (Figure  2D  too,  the  E-P  distribution  mimics 
the  mirror  image  of  P  w  ith  a  mean  of -1.40  mm/day  and  a 
standard  deviation  of  2.43  mm/day.  The  flux  of  freshwater 


into  this  region  is  from  the  Indonesian  Throughflow  (ITF), 
Inter  Tropical  Convergence  Zone  (1TCZ)  rainfall  and 
freshwater  transported  from  the  BoB  [Joseph  and  Freeland , 
2005;  Sengupta  et  al. ,  2006].  E-P  is  highest  in  the  AS 
(Figure  2e)  with  a  mean  of  3.47  mm/day  and  a  standard 
deviation  of  2.17  mm/day.  This  reflects  the  minimal  influ¬ 
ence  of  freshwater  in  the  distribution  of  SSS  in  this  basin. 

3.2.  Surface  Flow  Velocity  Structure 

[23]  In  this  section,  we  describe  the  structure  of  the 
HYCOM  simulated  surface  flow  velocities.  The  zonal 
component  shows  significant  seasonal  variability  in  the 
study  regions  except  in  the  EQ  region,  where  it  is  minimal 
(Figures  3a-3d).  Surface  flow  is  eastward  during  the  SW' 
monsoon  at  which  time  the  maximum  velocities  are  also 
observed.  Surfaee  flow  is  westward  during  the  NE  monsoon 
in  all  the  study  regions  cxeept  the  EQ  region,  where  it  is 
westward  during  the  SW  monsoon  months  and  becomes 
eastward  (decrease  of  negative  westward  velocity)  by  April/ 
May  under  the  influence  of  equatorial  jets  (Figure  3d).  The 
meridional  component  of  the  flow  shows  significant  seasonal 
variation  only  in  the  SL  region  (Figure  3e).  Overall,  the 
strongest  velocities  are  seen  in  the  SL  region  (Figure  3c) 
emphasizing  its  role  in  the  exchange  of  salt  and  mass 
between  the  AS  and  BoB.  The  seasonal  reversing  of  these 
velocities  afTeets  the  currents  and  hence  the  transport  of 
water  and  salt.  It  also  contributes  to  the  seasonal  variation  in 
the  SSS  transport  that  is  observed  especially  in  the  northern 
Indian  Ocean. 

3.3.  Seasonal  Variation  of  SSS  Transport 

[24]  Results  for  meridional  and  zonal  surfaee  salinity 
transport  magnitudes  are  shown  in  Figure  4.  The  patterns  of 
magnitude  and  direction  of  the  surface  salt  transport  track 
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Figure  3.  Seasonal  variation  of  surface  velocity  flow  in  the  Arabian  Sea  (AS),  Bay  of  Bengal  (BoB),  Sri 
Lanka  (SL),  and  Equatorial  (EQ)  regions. 


that  of  the  surface  flow  velocity  (Figure  3).  Zonal  transports 
are  higher  and  show  larger  seasonal  variation  than  meridi¬ 
onal  transport  in  all  the  study  regions  (Figure  4).  In  the  AS 
and  BoB  regions  (Figures  4a  and  4b),  zonal  SSS  transport  is 
maximum  (0.75  *  106  kg  s”1)  and  eastward  during  the  peak 
of  the  SW  monsoon.  Transport  becomes  westward  during 
the  NE  monsoon  season  but  of  lesser  magnitude  (maximum 
of -0.2  x  106  kg  s“l)  than  during  the  SW  monsoon  season. 


The  westward  flow  in  the  BoB  during  winter  is  responsible 
for  the  transport  of  low-saline  water  via  the  North  Equatorial 
Current  into  the  AS.  In  the  SL  region  (Figure  4e),  maximum 
zonal  salt  transport  (0.5  x  106  kg  s“!)  occurs  at  the  peak 
of  the  monsoons  and  of  similar  magnitudes  but  reversed 
directions.  It  is  westward  and  eastward  respectively  during 
the  NE  and  SW  monsoons.  In  the  EQ  region  (Figure  4d), 
zonal  surfaec  salinity  transports  arc  almost  entirely  westward 


Figure  4.  Seasonal  variation  of  area-averaged  zonal  and  meridional  surface  salinity  transport  for  the  (a) 
Arabian  Sea  (AS),  (b)  Bay  of  Bengal  (BoB),  (e)  Sri  Lanka  (SL),  and  (d)  Equatorial  (EQ)  regions. 
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Figure  5.  Comparison  oflatitudinal  variation  of  net  meridional  depth-integrated  salt  transport  eomputed 
using  HYCOM  and  SODA. 


with  maximum  (-0.45  x  10*  kg  s-1)  occurring  in  September, 
at  the  end  of  the  SW  monsoon  season  while  the  minimum 
(0.05  x  10*  kg  s”1)  occurs  in  February,  at  the  end  of  the  NE 
monsoon  season.  Meridional  SSS  transport  is  northward 
in  the  AS  during  the  NE  and  SW  monsoon  seasons  with 
maximum  (~0.1  *  10*  kg  s-1)  observed  at  the  peak  of  the 
monsoons.  It  is,  however,  southward  during  the  pre-SW 
monsoon,  where  it  reaches  -0.1  x  106  kg  s_1  in  May.  In  the 
BoB,  meridional  transport  is  mostly  southward  during  the 
SW  monsoon  season  and  at  the  start  of  the  NE  monsoon 
season  but  reverses  to  southward  at  the  peak  of  the  later 
season.  The  largest  meridional  SSS  transport  occurs  in  the 
SL  region,  where  it  is  northward  during  the  NE  monsoon 
season  with  a  peak  of  about  0.1  x  106  kg  sTl.  Dunng  the  SW 


monsoon  season,  the  transport  reverses,  becoming  south¬ 
ward  and  peaking  at  about  -0.2  x  10*  kg  s  1  toward  the  end 
of  that  season.  In  the  EQ  region,  meridional  SSS  transport 
fairly  tracks  that  of  the  zonal  transport  and  is  mainly  south¬ 
ward  with  northward  transport  observed  only  at  the  end  of 
the  NE  monsoon  season. 

3.4.  Depth-Integrated  Meridional  Salt  Transport 

[25]  Depth -integrated  meridional  salt  transport  computed 
using  HYCOM  arc  compared  to  that  using  SODA  (a  data 
assimilated  product)  with  the  results  (Figure  5)  showing 
they  track  each  other  quite  reasonably  in  sign  and  magni¬ 
tude.  The  annual  mean  estimates  of  HYCOM  net  meridional 
salt  transport  in  the  Indian  Ocean  (Figure  6)  indicate  a 
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20  40  60  80  100  no 

Longitude 


Figure  6.  Annual  mean  depth -integrated  net  meridional  salt  transport  (x  10*  kg/s)  for  eaeh  5°  latitudinal 
belt  in  the  Indian  Ocean. 
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Figure  7.  Latitudinal  variation  of  depth-integrated  meridional  salt  and  volume  transport  in  the  Indian 
Ocean  computed  from  4  year  mean  HYCOM  simulations. 


southward  transport  at  all  latitudes  except  20°N,  where  it 
is  396  x  1()6  kg/s.  It  generally  increases  southward  with 
the  highest  transports  (-463  *  106  kg/s  to  -531  *  106  kg/s) 
occurring  in  the  south  (10°S-35°S)  and  a  maximum  of -552 
x  106  kg/s  at  30°S.  The  northern  salt  transport  at  20°N  may 
be  due  to  the  shallowness  of  the  bottom  topography  and  all 
the  layers  would  not  be  considered  as  that  at  for  example 
20°S  or  30°S,  where  the  ocean  is  deep.  This  northward  salt 
transport  at  20°N  and  southward  salt  transport  at  15°N  in 
the  Arabian  Sea  would  lead  to  salt  divergence  between  15°N 
and  20°N  and  we  do  not  have  any  immediate  explanation 
for  this.  Probably,  this  is  an  important  feature  that  leads  to 
the  flow  of  high-salinity  water  masses  of  Persian  Gulf  origin 
and  Red  Sea  origin  at  subsurface  depths  and  also  compen¬ 
sates  the  net  salt  loss  due  to  predominant  evaporation  than 
precipitation  and  no  river  runoff.  The  southward  direction  of 
salt  transport  in  the  Indian  Ocean  is  consistent  w  ith  the  study 
of  Toole  and  Warren  [1993]  where  they  obtained  a  net  salt 
flux  of -236.8  x  106  kg/s  at  32°S. 

[26]  The  net  salt  and  volume  flux  track  each  other  quite 
closely  at  all  latitude  except  bctw'ccn  5°S  and  20°S,  where 
there  is  a  significant  difference  in  the  patterns  (Figure  7). 
This  is  as  a  result  of  large  volume  of  water  from  the  1TF. 
This  water  is,  however,  less  saline  and  therefore  does 
not  significantly  increase  the  magnitude  of  salinity  in  that 
region. 

[27]  Sectionally  integrated  salt  fluxes  showing  contribu¬ 
tions  from  each  layer  arc  shown  in  Figure  8.  In  the  northern 
Indian  Ocean,  at  20° N  (Figure  8a),  southward  transports 
arc  seen  only  in  deep  waters  with  the  rest  of  the  water 
column  showing  northward  flux.  The  net  flux  of  salt  at  this 
latitude  is  northward.  The  flow  is  southward  in  almost  all 
layers  at  10°N  (Figure  8b)  with  a  northward  transport 
between  layers  12-15.  Generally  for  the  Southern  Indian 
Ocean  (Figures  8d-8f),  surface  and  deep  water  layers  show' 
southward  transport  while  intermediate  and  bottom  water 
layers  show'  northward  transport.  Net  fluxes  for  the  entire 


sections,  however,  show  a  northward  salt  transport  of 
bottom  waters  and  a  southward  salt  transport  of  surface 
waters.  When  considered  as  a  whole,  the  Indian  Ocean 
shows  a  net  southward  transport  of  salt.  Plots  of  depth- 
integrated  volume  transport  (Figures  9a-9f)  show'  that  they 
track  that  of  the  salt- integrated  transport  closely  in  both 
magnitude  and  direction. 

[2k]  Freshwater  flux  mainly  directly  affects  distribution 
and  variability  of  surface  salinity  [Wi[ffels  et  a/.,  1992;  Rao 
and  Sivakumar ,  2003].  In  the  deeper  regions  of  the  Indian 
Ocean,  salt  distribution  is  influenced  by  water  transported 
into  it  from  other  oceans.  Deep  water  is  not  formed  in  the 
Indian  Ocean  but  rather  they  originate  from  the  Atlantic  and 
Southern  Oceans.  Localized  deep  water  is  only  formed  by 
mixing  [Mantyla  and  Reidy  1995].  As  a  result,  the  Indian 
Ocean  converts  cold  dense  water  to  warm  surface  waters, 
a  direct  opposite  to  what  happens  in  the  Atlantic  Ocean. 
Portions  of  the  southward  flowing  North  Atlantic  Deep 
Water  that  fails  to  be  incorporated  into  the  thcrmoclinc  of 
the  Atlantic  flows  eastward  with  the  Antarctic  Circumpolar 
Current  (ACC)  into  the  Indian  and  Pacific  Oceans  [ Gordon , 
1986;  Rintoul ,  1991].  Part  of  the  NADW  that  enters  the 
Pacific  Ocean  is  later  transported  into  the  Indian  Ocean 
via  the  ITF,  the  water  of  which  flows  westward  courtesy 
the  South  Equatorial  Current  (SEC),  all  the  way  to  the 
Mozambique  Channel  Along  this  transport  route,  large- 
scale  upwelling  returns  the  high-saline  deep  waters  to  the 
surface,  significantly  changing  the  salt  content  of  the  Indian 
Ocean  [ Stommel ,  1958;  Gordon ,  1986;  van  Aken  et  a/., 
2004]. 

[29]  The  ITF  plays  important  role  in  the  distribution  of 
salt  in  the  Indian  Ocean.  Its  strength  influences  the  salinity 
budget  as  far  west  as  the  Somali  Current  region  [ Godfrey , 
1996;  Gordon  et  a/.,  1999;  Potemra ,  2005].  In  this  study, 
we  obtained  an  ITF  volume  flux  of  14.4  Sv  into  the  Indian 
Ocean;  a  value  close  to  Piola  and  Gordon's  [1984]  14  Sv. 
Other  estimates  include  12  Sv  [ Godfrey ,  1989],  7  Sv  [Fu, 
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Figure  8.  Depth-integrated  meridional  salt  transport  in  eaeh  layer  integrated  along  latitudinal  sections  in 
the  Indian  Ocean  computed  from  4  year  mean  HYCOM  simulations.  Positive  transport  is  to  the  north. 


1986],  and  18.6  Sv  [Fieux  el  al.,  1994].  The  different  values  about  800  m  (Figure  10).  This  pattern  traeks  that  of  the  salt 

underscore  the  difficulty  in  estimating  the  1TF  as  a  result  flux  which  shows  a  net  flux  of  512.2  *  106  kg  s_l  into  the 

of  its  complex  topography  and  large  variation  in  tides  Indian  Oeean. 

[Potemra,  2005;  Valsala  and  lkeda ,  2007].  It  also  could  [30]  There  is  higher  transport  at  higher  latitudes  in  the 
imply  a  large  seasonal  and  interannual  variation  of  the  ITF  Indian  Oeean.  The  largest  seasonal  variation  (January  and 

[Godfrey,  1996].  From  the  HYCOM  estimates,  the  bulk  of  July  compared)  in  net  meridional  salt  transport  is  seen  eloser 

the  ITF  flux  occurs  in  the  top  18  model  layers  representing  to  the  equatorial  region;  between  10°S  and  20°N,  with  the 


I  ransport  (Sv) 

Figure  9.  Depth-integrated  meridional  volume  transport  (1  Sv  =  I06  m3  s'1)  in  each  layer  integrated 
along  latitudinal  sections  in  the  Indian  Oeean.  Positive  transport  is  to  the  north. 
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Figure  10.  Depth -integrated  zonal  (a)  salt  and  (b)  volume  transport  in  eaeh  layer  integrated  along  the 
ITF  region.  Positive  fluxes  are  to  the  east. 


largest  occurring  at  10°S  and  10°N  (Figure  11).  Currents 
southward  of  10°S  do  not  show  much  seasonal  variability 
as  compared  to  those  northward  of  10°S  [Schott  and 
McCreary y  2001].  At  the  southernmost  extent  of  the 
Indian  Ocean,  there  is  not  much  variability  in  the  salt  tlux. 
Between  15°S  and  35°S,  the  seasonal  values  are  almost  the 
same.  Between  15°S  and  15°N,  the  net  flux  is  southward 
and  July  salt  tlux  is  generally  lower  than  the  annual  mean 
salt  flux  while  January  salt  flux  is  generally  higher  than  the 


annual  mean.  From  20°N  to  25°N,  however,  the  net  flux  is 
northward  and  July  flux  is  higher  than  the  annual  mean  while 
the  January  flux  is  lower  than  the  annual  mean.  Flux  de¬ 
creases  as  latitude  decreases  from  the  south  to  the  north.  The 
northern  flux  is  weaker  as  the  effect  of  freshwater  flux 
“dilutes"  the  water.  For  the  entire  study  period  (2003-2006), 
the  monthly  mean  salt  transport  (Figure  12)  at  the  selected 
latitudes  south  of  the  equator  show  year-round  southward 
flux.  Northward  of  the  equator  and  especially  at  20°N, 


Figure  11.  Seasonal  variation  of  net  depth -integrated  meridional  salt  transport  in  the  Indian  Ocean 
computed  from  4  year  mean  January  and  July  simulations.  Positive  transport  is  to  the  north. 
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Figure  12.  Seasonal  variation  of  depth- integrated  meridional  salt  transport  in  each  layer  integrated  along 
latitudinal  sections  in  the  Indian  Ocean.  Positive  transport  is  to  the  north. 


northward  transports  are  observed  in  some  months.  At 
10°S,  20°S  and  30°S,  maximum  meridional  transport  (-30  x 
IO6  kgs  ')  occurs  in  June,  at  the  start  ofSW  monsoon.  During 
the  NE  monsoon  season,  the  transports  are  more  toward  the 
equator  while  being  more  toward  the  poles  during  the  SW 
monsoon  season. 

[31]  Our  estimates  of  meridional  Ekman  volume  and  salt 
transport  (Figure  13)  show  a  net  southward  transport  in  the 
Indian  Ocean.  This  trend  in  the  Indian  Ocean  has  been 
attributed  to  the  strong  summer  monsoon  southwesterly 


winds  [Levitus,  1988;  Chereskin  ct  al .,  2002].  The  volume 
transport  magnitude  and  direction  is  close  to  that  of  Levitus 
[1988],  who  showed  similar  magnitudes  and  directions. 
The  annual  mean  meridional  Ekman  volume  transport  is 
northward  (~0.4  Sv)  only  at  20°N.  Equatorward  of  20°N,  it 
becomes  increasingly  southward  between  reaching  a  maxi¬ 
mum  of  about  10  Sv  at  2°30rN.  No  data  is  presented 
between  2°30'N  and  2°30'S  as  a  result  of  the  absence  of 
Coriolis  near  the  equator.  In  the  southern  Indian  Ocean, 
meridional  Ekman  volume  transport  increases  poleward 
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Figure  13.  Annual  mean  latitudinal  variation  of  depth -integrated  Ekman  meridional  salt  and  volume 
transport  in  the  Indian  Ocean. 
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from  about  “0.5  Sv  until  it  reaches  a  maximum  of  about 
13  Sv  at  12°30'S.  Levitus  [1988]  reported  a  maximum  of 
-18  Sv  at  1 1°N.  Southward  of  12°30'S,  the  Ekman  merid¬ 
ional  volume  transport  is  still  southward  but  of  lesser 
magnitude.  The  meridional  Ekman  salt  transport  traeks  that 
of  the  volume  closely,  with  a  maximum  (25  *  10*  kg  s'1) 
occurring  at  12°S.  A  secondary  maximum  of  -24  *  106  kg 
s'1  occurs  at  2°30fN.  A  positive  Ekman  salt  transport  (-1  * 
10*  kg  s'1)  is  observed  at  2°30rN.  The  reduced  meridional 
Ekman  transports  observed  southward  of  12°30'S  can  be 
attributed  to  the  lesser  impact  of  monsoon  winds  in  that 
region. 

4,  Conclusions 

[3:]  In  this  paper,  4  year  HYCOM  simulations  are  used  to 
study  the  transport  of  salinity  in  the  surface  layers  of 
the  Indian  Ocean  and  to  estimate  the  flux  of  salt  in  the 
32  HYCOM  layers  using  depth-integrated  transport  methods. 
Results  show  the  temporal  and  spatial  variability  of  surface 
salinity  variability  terms  among  the  four  subregions  studied. 
Precipitation  is  highest  in  the  Bay  of  Bengal  while  showing 
a  sinusoidal  pattern  in  the  Sri  Lankan  region.  Although 
evaporation  is  high  in  the  Arabian  Sea,  seasonal  variation  is 
quite  minimal  in  the  study  regions  w  ith  annual  bimodal  peaks 
in  the  AS  coinciding  with  the  minima  and  maxima  respec¬ 
tively  ofE-P.  Zonal  SSS  transport  was  larger  than  meridional 
SSS  transport.  The  strongest  variability  is  seen  in  the  SL 
region  while  in  the  EQ,  SSS  transport  was  almost  entirely 
southward  and  westward  respectively.  Depth -integrated 
transport  studies  showed  a  northward  transport  of  volume 
and  salt  in  bottom  water  layers  and  a  southward  transport 
of  salt  in  surface  layers.  The  highest  net  salt  flux  of  —552  * 
10*  kg/s  is  observed  at  30°S,  a  region  influenced  by  the 
Cireumpolar  Deep  Water  and  the  Antarctic  Bottom  Water. 
Seasonal  variation  in  salt  flux  is  significantly  notieeable 
along  the  equatorial  region  where  1TCZ  influences  rain 
patterns  and  as  well  as  influence  by  the  1TF.  The  1TF 
transports  about  14.4  Sv  of  water  and  512.2  *  10*  kg  s“l  of 
salt  into  the  Indian  Ocean.  Overall,  there  is  a  net  southward 
flux  of  salt  in  almost  all  latitudes  of  the  Indian  Ocean  with 
patterns  reflecting  the  transport  of  water  masses  and  currents 
in  the  Oeean.  Analyses  of  meridional  Ekman  volume  and  salt 
transport  show  a  predominantly  southward  transport,  an 
indication  of  the  strong  influence  of  SW  monsoonal  winds. 
This  study  illustrates  that  the  HYCOM  output  is  useful  for 
calculating  Indian  Oeean  salinity  transport  and  should  be 
useful  in  computing  near-surfaee  salt  transport  from  the 
Soil  Moisture  and  Oeean  Salinity  (SMOS)  and  Aquarius 
satellite  missions  as  well. 
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